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Abstract 

We present a model for the electric potential profile across the membranes of neuronal cells. We considered the resting and 
action potential states, and analyzed the influence of fixed charges of the membrane on its electric potential, based on 
experimental values of membrane properties of the spinal ganglion neuron and the neuroblastoma cell. The spinal ganglion 
neuron represents a healthy neuron, and the neuroblastoma cell, which is tumorous, represents a pathological neuron. We 
numerically solved the non-linear Poisson-Boltzmann equation for the regions of the membrane model we have adopted, 
by considering the densities of charges dissolved in an electrolytic solution and fixed on both glycocalyx and cytoplasmic 
proteins. Our model predicts that there is a difference in the behavior of the electric potential profiles of the two types of 
cells, in response to changes in charge concentrations in the membrane. Our results also describe an insensitivity of the 
neuroblastoma cell membrane, as observed in some biological experiments. This electrical property may be responsible for 
the low pharmacological response of the neuroblastoma to certain chemotherapeutic treatments. 
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Introduction 

Electrostatic forces affect the passive and active transport of 
charged particles through biological membranes. The flow rate of 
ions through the membrane depends on the strength of the 
intramembranous electric field. These forces also affect the 
robustness of some ligands of the membrane [1]. In this work, 
we study the influence of surface electric charges on the stability of 
the cell membrane in the condition of equilibrium, by modehng 
the electric potential profile. The profile describes the behavior of 
the potential along the axis perpendicular to the cell membrane, 
from the outer bulk region to the inner cytoplasmic region [2-5]. 
We do not consider here dynamical phenomena in the structure of 
the membrane, and treat only the electrostatic situation, which 
occurs once the system has reached equilibrium. We refer the 
reader to studies such as [6,7] that treat dynamical, nonequilib- 
rium phenomena, like the molecular dynamics of ion channels 
associated with transmembrane ion transport, using the Poisson- 
Nersnt-Planck theory [6] and the Poisson-Boltzmann-Nernst- 
Planck model [7] . 

The electric potential on a cell surface is determined as the 
difference of potential between the membrane-solution interface 
and the bulk region [1]. It has been shown that the electrophoretic 
behavior of neuroblastoma cells provides information about their 
surface charges, in different phases of the cellular cycle [8—10]. 
These experiments show that membrane anionic groups are 
mainly responsible for the surface charges of murine neuroblas- 



toma cells [10]. It is known that neuroblastoma cells, like all other 
cancerous cells, multiply quickly. Alterations of the dynamics of 
cellular multiplication compromise the synthesis and structure of 
components of the membrane, with possible degradation of these 
components, promoting deformations of the structure and 
composition of the plasma membrane [11]. 

We show a detailed and revised description of the model more 
briefly presented by Cortez and collaborators in [3-5] , which was 
originally used to simulate the squid giant axon. This model is 
based on the statistical mechanical theory of electrolyte solutions 
and electric double layers [12-15]. We then present a study that 
applies this model in a novel way to the neurons of mammals 
(mice) [16,17], in order to investigate the alterations of the electric 
potential and therefore, the capability of transmitting electric 
signals in the membrane of cancerous neurons. Here, the spinal 
ganglion neuron denotes a healthy neuron, and the neuroblastoma 
cell represents a tumorous neuron. With simulations of this model, 
we compare the effects of charges fixed onto the inner surface of 
the membrane and those associated with cytoplasmic proteins, on 
the electric potential on the surfaces of the membranes of both 
types of cells, considering both natural states of neurons, the 
resting and the action potential (AP) states. The AP state refers to 
the state of the neuron in which it has been stimulated enough, so 
that its physico-chemical conditions are such that the transmem- 
brane potential reaches the maximum value of the AP. The 
temporal evolution of the transmembrane potential was not 
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considered. We also calculated the potential profile across the 
membrane, including data from electrophoretic experiments in 
our model. 

Methods 

Cortez and collaborators have proposed in [4] an axon 
membrane model, to study how charges fixed onto the inner 
surface of the membrane and those associated with cytoplasmic 
proteins influence the electric potential of the squid axon 
membrane. In their work, the effects of divalent ions were 
included, with a numerical solution of the model equations. In the 
present study, we apply this model to healthy and cancerous 
mammalian neurons, to understand electrical characteristics of the 
membranes of these cells. We present the formal derivation of the 
complete model here, since it was not shown in [3,4] and because 
we have found mistakes in some of the equations presented in [4] . 
We also describe, in Section Surfcwe Potentials, a method for 
calculating electric potentials on the interfaces of the neuronal 
membrane, which was not discussed in previous work. We thus 
present here a more detailed and revised description of the theory 
and mathematical model of the electric potential across neuronal 
membranes in equilibrium, which was originally discussed in [4]. 

In the neuronal membrane model we have adopted, shown in 
Fig. (1), four different regions are presented: extracellular, 
glycocalyx, bilayer and cytoplasm. The bilayer thickness is h and 
the width of the glycocalyx is hg. Surface potentials are represented 
as (^5^^ for the potential on the surface S^g, between the 
extracellular and glycocalyx regions, i^^^^^ is the potential on the 
surface Sgb, between the glycocalyx and the bilayer, and is the 
potential on the surface She, between the bilayer and cytoplasm. 



We denote by and 
eletlrolytic extracellular phase, 
cytoplasmic region, respectively. 



the potentials at — oo, in the 
and at + oo, in the bulk 



The Electric Potential in the Membrane Regions 

To determine the potential profile across the membrane, we 
first consider the Poisson equation [3,4,15], 



-4jrp,(z) 



(1) 



for / e {ext,g,b,c}. 



where <l>i(x,)',z) is the electric potential in any region /; i = ext for 
the outer electrolytic region; i = g for the glycocalyx; i = b for the 
bilayer; and i = c for the cytoplasm. The volumetric charge density 
due to the electrolytes in solution of area ; is and £,■ is the 

dielectric constant in region There is no charge in the bilayer 
{i = h), due to its hydrophobic property, and thus Pi,iz) = Q. 

We consider the further boundary conditions to calculate the 
constants of integration: 



1 . When z tends to an extreme 
z— » + 00 (in the cytoplasm). 



'alue, z— > — oo (in region ext) or 
the electric potential tends to 



limiting values represented by 



and 



respectively. 



Ionic concentrations assume limiting values in the bulk regions, 
lim '?y,«((^) = ';y,-oo and lim »j (z) = >7 _|_„, respectively, 

where y represents an ion, such as Na"*", K"'", or CI~. 

2. When z = —hg — h/2 and z= +h/2, we consider the continu- 
ity of the electric potentials, (j>g^^, (j) and (j) g^^^ , respectively. 

3. There is a discontinuity of the electric field vector on the 
surfaces between regions. 



The Effect of Fixed Charges. When we include the effect of 
fixed charges in the model, Poisson Eq. (1) becomes 



V^iix,y,z)- 



(2) 



where py, is the density of charges fixed onto proteins of area / 
(Pfext ~ 0) for the outer electrolytic region, pjg for the glycocalyx, 
and pj^ for the cytoplasm). 

The volumetric charge density p,(z) is the sum of the charge 
densities of positive and negative ions in the solution [3] 



(3) 



where y+ represents a positive ion, and y , a negative ion. The 
molar density for an ion y (ionic concentration) in region i is t]y ,(z), 
and Vy is the valency of ion y. For example, V^^2+ =2 and 
Kc/- = — 1. The absolute value of the electron charge is e. Due to 
the electroneutrality condition, we can write 



J Pi{z)dz = 0, 

within the boundaries of ('ach region. 

In our model, we suppose that the surfaces are infinite in the x 
and y directions, perpendicular to z, and that the distribution of 
charges in these directions is homogeneous. Considering this, and 
substituting Eq. (3) in Eq. (2), we obtain 



dz^ 



where 



471 



Anpfl 



(5) 



To determine t]yj(z), we use the equation for the electrochemical 
potential, due to an ionic solute in a diluted solution [3,4] 



/i,,,(z) = nl^P,T)+KTln{ny,i{2)) + eVyUz), 



(6) 



where K is Boltzmann's constant, T is the temp(;rature, i^,f(P,T) 
is the standard chemical potential, dependent on pressure and 
temperature, KTln{r]y i{z)) is a term that expresses the influence of 
the ionic concentration r]yi{z), and eVy(j>i(z) is the contribution of 
the electric potential. 

Applying the V operator in Eq. (6) and again considering the 
homogeneous distribution of charges in the directions perpendic- 
ular to z, we verify that 



dHy^^iz) _ dnl,(P,T) 



d^i(z) 



dz dz ^KT-HnyA^)HeVy^. (7) 



Considering that there is a condition of Boltzmann equilibrium in 
the aqueous environments adjacent to the bilayer. 



%^=Oand 
dz 



= 0, and Eq. (7) becomes 

dz 
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Figure 1. Model for a neuronal membrane. Different regions are presented, with the corresponding symbols for the potentials on the surfaces 
dividing regions. Symbols are explained in the text. Minus signs illustrate negative fixed charges on proteins. 
doi:1 0.1 371 /journal.pone.00961 94.g001 



0 = KT^HU.)) + ev/-^. 



(8) 



The integration of Eq. (8) from z in one of the three regions, 
extracellular, glycocalyx and cytoplasm, to a limiting boundary 
region, for which we have experimentally measured quantities, 
considering the electrolytes distributed over the adjacency of the 
bUayer gives 



where i^,- / and t]y jj are limiting values of the electric potential and 
the ionic concentration of y, respectively, in region ;'. The solution 
of Eq. (9) results in 



ny,i(z) = riyjj exp 



where 



A0,(z) = ^,-(2)-f-,/- 
The molar density for a positive y ion is thus given by 



and for a negative y ion 



ny-,iiz) = ny-,i.i exp 







(10) 



(11) 



(12) 



(13) 



Zy = \V;\ 



(14) 



Equations (12) and ( 1 3) are the Boltzmann distribution of charges 
due to the presence of positive and negative y ions [15], 
respectively, in the phases adjacent to the bilayer. Substituting 
Eqs. (12) and (13) in Eq. (4), we obtain 



4n 



(15) 



'^eZy-riy- jj exp 



) 



In the bulk regions, we can consider the electroneutrality condition 

^ eZ,,+ 1?.,+ i i = J2 eZy- riy-jj, (16) 

1'- 

and, in a first approximation, we assume a symmetric electrolyte to 
simplify our calculations, so that 



1y+.iJ = 1y-,iJ = 1yJJ- 



(17) 



We have taken the boundary values from experimental measure- 
ments in the bulk regions and on surface Seg, so that (l>ext,l = - oo > 
and <l>g i = (l>s„„- For the ionic concentrations. 



where 



'?•,.,<.«,/ = '?■/,- 00 > 1y,cj = ^y, + ao, and r]ygi = r]yg(-hg-h/l). 
Throughout, we denote by ( 

Sij between regions / and / 



Throughout, we denote by ^g.. the electric potential on surface 
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We can now use Eq. (17) to rewrite Eq. (15) as 



X exp 



If we consider that 



\ KT J \ KT J 



■ exp 



(18) 



^^^^-exp^^^=-2sinhf^^ 



exp — exp 

we can write Eq. (18) as 



V KT 



(19) 



eZy^iiz) 



KT 



+ai. (20) 



Our model only considers mono (Zi = 1) and divalent (Z2 = 2) ions 
[4]. We thus limit the y sum to 



dz^ 



87re;j,,-, 



+ 



(21) 



To simplily further calculations, we denote 



ti KT' 



(22) 



and Eq. (21) may be expressed as 



=Ai sinh(M0,(z)) +5; sinh(2M<^,.(z)) (23) 
Considering that 

sinh(2^A(^,.(z)) = 2 sinh(j8A<^i(2)) cosh(jSA0,.(z)) , (24) 
Eq. (23) can be rewritten as 



d^i(z) 



--A.smh{tSA,l>,(z)) 



dz^ -—-^r-r,.-., ^25) 

+ 2Bi sinh(j8A0,.(z)) cosh(jSA0,.(z)) +0,-. 



Multiplying Eq. (25) by 



2c^^, = 2^Jz, 
dz 



(26) 



and integrating, we have 



J^,(z)d<^^(z) ^_ 
dz^ dz 



=|2^,sinh(jSA(^,(z))(/,^,- 
+1 4Bi sinh()SA(^,(z)) cosh(M<^,(z))«^'^i (27) 
+ I luid^i. 



whose solution is 



dz 



_ 2Ai cosh(;6A(^,.(z)) 25; cosh2(Mi^,(z)) 
^ P ^ ^ (28) 

+ 2ai(^,.(z)+fF,-, 



where Wi is a constant of integration for region /. 

To further simplify the symbolic representation of the equations, 
considering Eqs. (5) and (22), we denote 



Qmi=2Ai, Qd.=2Bi, and, 
gi = 2ai. 



(29) 



where 2^,- and Qj. express the influence of boundary value 
concentrations of monovalent and divalent ions, respectively, in 
each region Eq. (28) is thus expressed as 



dm 

dz 



2 e„. cosh(M0,(^ ^ erf.cosh'(M^,(z)) 



(30) 



+gi(l>i(z)+Wi. 
We can rewrite Eq. (30) as 



where 



'^=Vm+g4Az)+m, 



Qm, cosh(jSA^,(z)) ^ Qdi cosh^ (iSA^,(z)) 
m= -^—^ + ^ J . 



(31) 



(32) 



Equation (31) is the Poisson-Boltzmann equation for the electric 
potential in any region /. 

Solution of die Poisson-Boltzmann Equation for die 
Extracellular Region. {i = exf) 

In the extracellular region, the effect of fixed charges is 
negligible (p/ex( = 0 and gex; = 0), and the solution of Eq. (31) 
therefore only considers the electrolytic charges. Moreover, the 
electric potential in z = — 00 is constant and we can write 



dz 



(z=-co) 



= 0. 



(33) 



In order to use Eq. (31) to calculate the potential in the 
extracellular region, we must find the value of W^y^t. We thus 
consider an imaginary surface 5_co, perpendicular to the z — axis, 
at zi<fi{ — hg — hl2), where Pfext = ^- We then integrate Eq. (31) 
from another position Z2 < Zi to Zi . Since both Z\ and Z2 are in the 
-00 region, <l>-oo=<l>exii^l) = <l>exii^2) = (l>s-ao' and 
A^gj.,(z2) = ^(,j.,(z2) — =0- can then substitute Eq. (33) 
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and A0g^,(z2) = O in Eq. (31), to calculate Wext as 



The constant of integration Wc is 



0 = 2^ cosii(O) + cosh2(0) + Wex,, (34) 



and 



imext ^aext 



(35) 



We then substitute Eq. (35) in Eq. (31) to obtain the differential 
equation for (j)^,^i(z), for any position z<{ — hg — h/2), in the 
extracellular region 



= V ^extKZ) + Oext(?), 



(36) 



where 



^ext{^) — 



e^,„[cosh(M0..X^))-l] 
i5 



(37) 



and, ^exti^) = <t>ext(^) -<t>- 
Equation (36) can be simplified to 



d't>exti^) / TTT 7T 

— = VXexti^) + (Oextiz), 

where 



(38) 



Xexti^) = 2 

and, 

COext(z) = 



1^ 



(39) 



Solution of the Poisson-Boltzmann Equation for the 
Cytoplasmic Region. (/ = c)Because, in the cytoplasmic re- 
gion, the potential in z = + oo is also constant. 



dz |(^=+")' 



= 0. 



(40) 



As we did for the extracellular region, to calculate the potential in 
the cytoplasmic region using Eq. (31), we must first find the value 
of Wc. We thus consider an imaginary surface S+co, perpendic- 
ular to the z — axis, at z^»(h/2). We then integrate Eq. (31) from 
another position Z4>Z3 to Z3. As both Z3 and Z4 are in the -I- go 
region, 0 ^ ^ = ^^(zj) = 0^(z4) = (^5^ ^ , and 

A(l>c(z4) = (l>c(z4)-(l>s^^=0. We thus substitute Eq. (40) and 
A^^(z4) = 0 in Eq. (31), to obtain fVc as 

0= %cosh(0) + %cosh2(0) +gc<l>+^ + W,. (41) 



j8 



- Qmc - Qdc 



-gc(l>+c 



(42) 



We then substitute Eq. (42) in Eq. (31) to obtain the differential 
equation for <j>ciz), for any z>{h/2), in the cytoplasmic region 



= V2,(z) + (T,(z) + g,A^,(z), 



(43) 



where 



2„,[cosh(M<^,(z))-l] 



grf,[cOSh"(M^c(^))-l] 



(44) 



and , A(l>ciz) = <l>ciz)-<l>+^. 
Equation (43) can be simplified to 



^ = VZc(2) + «c(z)+gcA^,(z), 

where 



(45) 



Qmc smh^'^A4>c(z) 



Wciz)-- 



p 

Qjc sinh^iPAU^)) 
P 



, and. 



(46) 



Solution of the Poisson-Boltzmann Equation for the 

Glycocalyx Region. (/ = g^)We consider that the potential <j>g^ 

on the surface Sg^ satisfies 



d(l>g(z) 
dz 



-A/2) 



-glScg 



(47) 



where z is the unit vector in the z direction, and Eir I is the 

u 

electric field in region k, Eki^ij), at Zy, the position of surface Sy 
between regions / and j, e.g., Eg\g^^=Eg{ — hg — h/7) is the 
electric field in the glycocal^'x region at the position of surface Seg. 

We have previously stated that the boundary condition for ionic 
concentrations in the glycocalyx is defined as 
riy^l = rjyg{ — hg—h/2). In order to obtain r]yg{ — hg — h/2) from 
the experimentally measured »7y__oo! we would need to apply Eqs. 
(12) and (13). This would result in four values oi rjy g( — hg — h/2) 
(for positive and negative, monovalent and divalent ions). In this 
case, the symmetric electrolytes assumption (Eq. (17)) would not 
hold, and the mathematical formahsm leading to Eq. (31) would 
not apply. Nevertheless, in Eq. (15), the difference in the 
contributions of the terms involving rj^ g i for the different ions is 
small compared to the value of the a,- term, which is a few order of 
magnitude larger. We thus assume that rjyg{ — hg — h/7) = rjy _^, 
for monovalent and divalent ions. 

In order to use Eq. (31) to calculate the potential in the 
glycocalyx region, we need to find the value of Wg. We thus solve 
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Eq. (31) at z= —hg — h/2, on the surface Seg, and take 
^gi-hg-h/^) = h,^, and 
A(j>g( — hg — h/2) = (l>g( — hg — h/2) — <j>g^^=0. We can then substi- 
tute Eq. (47) and K^g{ — hg-h/2) = Q in Eq. (31), to calculate Wg 
as 

{E,\seg)^= ^COSh(0)+ ^COSh^ {(i)+gg^Seg + ^g. (48) 

The constant of integration Wg, for the glycocalyx region, is 
therefore 



Surface Potentials 

In order to solve the differential Eqs. (38), (45) and (52) for the 
extracellular, cytoplasmic and glycocalyx regions of the neuronal 
membrane, respectively, we must know the values of the surface 
potentials ^Seg' 't'Sgt, 'I' Sic' Although membrane surface 

potentials in cells cannot be measured experimentally, it is possible 
to obtain analytical predictions for the values of ij)^^^ cind(j)g^, 
from basic electrostatic relations, as we wiU now show. 

Considering the surface densities of electric charges. Gauss' law, 
and the discontinuity of the electric field vector on the surfaces Seg, 
Sgb and She, we obtain 



2 (Q mg + Qdg 



(49) 



(^gEg\se.-'^ex,Eext\se, =4nQSeg^^ 



(56) 



We then substitute Eq. (49) in Eq. (31) to obtain the differential 
equation for <^g(z), for any position { — hg — h/2)<z<{—h/2), in 
the glycocalyx region 



d(t>g(z) 
dz 



- ^kg{z) + ag{z) +ggA,^g{z)+Eg\l^^, (50) 



<^bEb\s„.-<^gEg\s=^-!tQsgb2, and. 



(57) 



where 



e^^[cosh(M'^^(z))-i] 

^g(?)= 

2rf,[cosh2(M^^(z))-l] 



and , A(^^(z) = <^j,(z) - ^^^^ . 
Equation (50) can be further simplified to 



dz 



S |2 



where 



Z^(z) = 2 



em^sinh'(^A«i^(z) 



cogiz)-- 



j6 

e^^sinh^(;gA^g(z)) 



, and. 



(51) 



(52) 



(53) 



Solution of the Poisson-Boltzmann Equation for the 
Bilayer. {i = b) As mentioned earlier, because the bilayer is 
highly hydrophobic, pj(z) = 0, and therefore Eq. (2) assumes the 
form 



respectively, where Ei,\ 



^bc 



(58) 



--Eh (see Eq. (55)). In the 
above, Qseg> QSgi, and Qs^^ stand for the charge density on the 
Seg, Sgi, and Sbc surfaces, respectively. 
As 



(59) 



in order to determine the discontinuity of the electric field vector 
on the surfaces 5^^, Sgb and Sbc, we substitute Eq. (59) in Eqs. (56), 
(57) and (58), and obtain 



4nQseg-eeJ-^ 



■^glSei 



dz \(z=-hg-hl2) ^ 

■^5 



(60) 



47t2s 



d<S,g(z) . 
^« dz \(z=-h/2) 



z, and. 



(61) 



V2|^fe(x,3;,z) = 0, 



(54) 



and its solution is a family of linear functions. The electric field 
within the bilayer (see Eq. (59)), Eb, can thus be expressed as 



'I'Sgb 't's^e ^ 

Eb=^—, z. 



(55) 



-^^QSbc 



d<^e(^) I 
dz l(z=A/2), 



(62) 



As we have previously obtained the expression that determines the 
electric field within the bilayer, Eb, we substitute Eq. (55) into Eqs. 
(61) and (62), and obtain 



PLOS ONE I www.plosone.org 



6 



May 2014 I Volume 9 | Issue 5 | e96194 



Modeling the Electric Potential across Membranes 



Table 1. Values of simulation parameters for both the spinal ganglion neuron and the neuroblastoma cell. 



Parameter 


Symbol 


Value 


Value In CGS 


References 


Dielectric constant in region h 


£6 


2 


2 


[3,4] 


Dielectric constant in region / {i = ext,g,c) 




81 


81 


[4] 


Glycocalyx width 


ih 


2.5 nm 


2.50 X 10"' cm 


[4,9,16] 


Bilayer thicl<ness 


h 


7.5 nm 


7.50 X IQ-' cm 


[3,4,9,16] 


Concentration of monovalent ions in bulk extracellular region 




0.154 M 


9.27 X 10" cm--' 


[16,17,29,30] 


Concentration of monovalent ions on S^-g 


•Ua, 


0.154 M 


9.27 X 10" cm"-' 


[16,17,29,30] 


Concentration of divalent ions in bulk extracellular region 


111. -a 


0.002 M 


1.20 X 10'* cm~' 


[16,17,29,30] 


Concentration of divalent ions on S^g 


'h.S„ 


0.002 M 


1.20 X 10'* cm-' 


[16,17,29,30] 


Concentration of monovalent ions in bulk cytoplasmic region 


^l, + oo 


0.154 M 


9.27 X 10" cm"' 


[29-31] 


Concentration of divalent ions in bulk cytoplasmic region 


'/2, + oo 


0.0004 M 


2.41 X 10" cm"' 


[29,30] 


Potential in — oo, the extracellular region 




0 mV 


0 statV 


[3,4] 


Temperature 


T 


310 K 


310 K 


[4] 


Boltzmann's constant 


K 


1.38x10- 


J/K1.38X 10"" erg/K 


[32] 


Absolute value of electron charge 


e 


1.60x 10" 


" C 4.80 X 10"'° state 


[32] 


Viscosity of region i 


>'( 


0.1 Pas 


1 P 


[3] 



'CGS is the centimeter-gram-second system of units. 
doi:l 0.1 371 /journal.pone.00961 94.t001 



h ( dSJz) 

= ^ V - 7, r - % ^ -,,2, ) , (63) 



^*=^^*.-!(^'^2.,.+..^l,../2,), (64) 
respectively. 

Substituting Eqs. (52) and (60) in Eq. (63), we determine the 
expression to calculate the surface potential i^^^^ 



where 



h.=K ^ + t^' 



(65) 



K = yJh/2) + oj,(h/2)+gMAh/2l 



(68) 



and here, A0^(A/2) = <^5^^^ -(^4 + ^. 

The electric potential on the surface S^g, (l>s^g> is determined 
from electrophoretic experiments [9,16,17] and the Helmholtz- 
Smoluchowski equation [3] 



4nvi ' 



(69) 



where ji is the electrophoretic mobility, = ^Seg~^ 



is the zeta 

potential, and v,- is the viscosity of region /. As, in our model, we 
define the potential ^_^=Q, the surface potential (^e is 



47IV,/< 



(70) 



^ = lgi- + cojf ( - h/1) + gg^<llg{ - h/2) + 



•eg ^ext ^, \{z=-hg-hlT) 



(66) 



and here, A(j)g( — h/2) = (j>g ^—(j>g^^^. In the same manner, 

substituting Eq. (45) in Eq. (64), we determine the expression to 
calculate the surface potential (j). 



'gb 



-VK, 



(67) 



where 



Model Calculations for Spinal Ganglion Neurons and 
Neuroblastoma Cells 

We have used data obtained from experimental observations 
[9,16,17] for values of parameters, in order to solve the first order 
ordinary differential equations, obtained from the Poisson- 
Boltzmann Eq. (31), for the different regions of the membrane. 
Some experimental values were obtained from electrophoretic 
experiments. Since each kind of cell presents a specific electro- 
phoretic mobility, the values of some parameters are different for 
the spinal ganglion neuron and the neuroblastoma cell, in our 
calculations. Tables (1) and (2) show all experimental values of the 
parameters used to solve the equations for the ganglion and the 
neuroblastoma. The difference i^ + oq— <^_oo is called the trans- 
membrane potential and is denoted as in the resting state, and 
in the AP state. We have defined (^_^ =0 in our calculations. 
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so that ^jj = 0 + oo in the resting state, and <j>A= 4 + 00 ™ AP 
state. We have thus examined the influence of parameters that 
represent electrical properties of the membrane, on the resting and 
AP states, and analyzed the differences between the healthy spinal 
ganglion neuron and the neuroblastoma cell. 

We implemented an approximate heuristic for finding roots of 
functions, to calculate (fig^^^ and (j>g^^ from Eqs. (65) and (67), which 
is specified in the next subsection. As mentioned earlier, the 
potential <l>g^^ was calculated with Eq. (70), from data obtained 
from electrophoretic experiments. 

As we included the density of charges fixed onto proteins within 
the membrane regions in the Poisson Eq. (1), we obtained a non- 
linear Poisson-Boltzmann Eq. (31), whose analytical solution has 
not been found. We therefore calculated values of the potential 
profiles with Eqs. (38), (45) and (52) numerically, using the Runge- 
Kutta method. The model simulation code is available on GitHub 
at https:/ /github.com/pintotm/PLoSOne2014. 

Roots of the System of Non-linear Equations. Because 
there is no experimental method to directly measure values of the 
surface potentials 0^^^ and <t>s^! we use Eqs. (65) and (67) that form 
a system of non-linear equations with two variables, ^ and , 
to determine these values. Some research work [3,4,9,16,17] 
indicates that the values of these potentials in real cells are in a 
limited region of the ^5.^ ^05^ plane. This means that we are 
looking for one of the roots of Eqs. (65) and (67), in a known 
region. 

We can write the system of Eqs. (65) and (67) as 

^^,>, = g{Kc) ^^^^ 

Our simple method for finding the roots of this .system consists in 
the steps described in Algorithm_Roots. At the resting state, 
min = — 250 mV and max = — 50 mV, while, at the AP state, 
min =10 mV and max = 100 mV. 

Algorithm. Algorithm Roots 

Step 0: 3^^^ = 10"* mV 

Step 1: Choose an initial value for <j)g^^ 

in min < ^Stc max 

Step 2: =20 mV 

Step 3: WhUe > do 
begin 

If l<^Sfc-«!'^*s»J < -5^ then 
else 

return to Step 1 
end if 
end While 

Algorithm for determining (j>g^^ and i^Sj j for the spinal 
ganglion neuron and the neuroblastoma cell. With the 
4si^ value found with this procedure, we use Eq. (67) to obtain the 
corresponding ^Sg, value. The surface potentials and ^Sic ^'^^ 
the spinal ganglion neuron and the neuroblastoma cell, whose 
values are shown in the following section, were obtained with 
Algorithm_Roots. 



Results 

The hydrophflic heads of the phospholipids that form the 
surfaces of the bUayer are negatively charged or polarized, and 
attribute a fixed charge density to surfaces Sgi, and She- In the 

situation of electrostatic equilibrium, which we are analyzing, both 
bilayer surfaces, Sgi, and She, are surrounded by a "diffuse electric 
layer" formed by the motion office ions in the fluid ionic solution, 
under the influence of electric attraction and thermal motion [12- 
15]. 

It is known that the inner surface charge density (|2s,,J) of the 
membrane is significantly higher than the outer surface density, 
due to the presence of negatively charged heads of phospholipids 
on the inner surface (phosphatidylserine), while on the outer 
surface the presence of neutral phospholipids dominates [18]. 
Moreover, the net charges fixed onto cytoplasmic proteins (Ipfd) 
are considered to be higher than the fixed net charges distributed 
in the glycocalyx region [19-21]. 

With our mathematical model, we first investigate the effect of 
Qs,, and pj-^ on the electric potential on the surfaces of the 
neuronal membranes of the spinal ganglion and the neuroblasto- 
ma. However, there is little information in the literature, regarding 
experimentally obtained quantities related to electric charges fixed 
within biological membranes. This is mainly due to the difficulties 
involved in obtaining the experimental measurements of these 
quantities. We therefore examine the behavior of the potentials 
^5^^ and ^5^^ for a range of values of the ratios QsbdQsee and 
Pfcl Pfgy given known experimental values of Qs^^ and Pfg, i.e. Qs,, 
and Pf^ are multiples of Qs^^ and pjg, respectively. We note that 
both glycocalyx and cytoplasm and their surfaces are negatively 
charged, so that 2s,j<0, gs„ <0, pjg<Q and Pf,<Q. 

Figure (2) shows the behavior oi^Sgt and ^Si, as we increase the 
negative charge on She, i.e. as we decrease We notice that for 
both the resting and AP states, a decrease of 2sfc has almost no 
effect on the surface potentials of both neuronal membranes. 
These variations in Qs^^. only determine a small gradual decay of 
4 Sic during the resting state of the cells. The resting and AP states 
were specified by boundary conditions, i.e. specific parametric 
values applied to the model. Values for ^ s^, and (j) are different 
between these types of ceUs, due to their specific membrane 
properties. We also observe that 0^^^ remains constant at — 34.97 
mV and —25.17 mV, respectively for the spinal ganglion and 
neuroblastoma cells, when they switch from the resting to AP 
states, and vice versa. Moreover, during the AP state, (j) assumes 
values near the transmembrane potentials (Fig. (2B)). 

We also examine the electric potential on the surfaces of the 
membranes of the healthy and the cancerous cells, in response to 
variations in the density of charges fixed onto proteins of the 
cytoplasm. Figure (3) presents the resulting (l>g^^ and (^^^^ as we 
increase this density of negative charges in the cytoplasm, i.e. as we 
decrease pf^. For both the resting and AP states, (^^^^ remains 
constant when values of pj-^. vary. However, a decrease in pj-^ 
causes a substantial fall of (^j^^, for both types of cells, at resting 
and AP states. At the AP state, Sg^^ presents a quick drop when 
Pfi./ Pfg<20 and tends to an asymptotic value, for decreasing 
values of pj-^, for both types of cells. 

Besides investigating the effect of fixed charges on the potential 
on the surfaces of these membranes, we study how the electric 
potential profile changes across the membranes of spinal ganglion 
neurons and neuroblastoma cells, for the resting and AP states. 
Although the values of Qgf^^ and p^-^ are not known, 1 2^^^ | and \pj-^\ 
are much larger than the corresponding charges in the outer 
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regions. We thus chose fixed values of Qs^^ = ^QQs^^ and 
Pfc = 20pjg (which are the same fixed values in Figs. (3) and (2), 
respectively) to calculate the potential profile in Fig. (4), for the 
resting state and in Fig. (5), during the AP state. For both natural 
states of these cells, we verify an accentuated decrease of the 
potential along the z axis, from the extracellular region to the 
surface of the glycocalyx. This decay is slightly more substantial for 
the neuroblastoma than for the spinal ganglion neuron, although 
the shapes of both curves are very similar. 

In the spinal ganglion neuron, the potential maintains a 
substantial drop across the glycocalyx. Interestingly, this phenom- 
enon does not occur for the cancerous cells, whose electric 
potential remains almost unchanged in this region. We observe 
linear variations of the potential across the bUayer of these cells, for 
both resting and AP states. During rest, the electric potential 
assumes much lower values on She than the asymptotic value in 
the bulk cytoplasmic region, while during the AP state, these 
values are similar. Before reaching the transmembrane potential 
value, the intracellular potential exponentially increases from the 
inner surface of the membrane to the bulk cytoplasmic region, 
during the resting condition, while during the AP state, we see no 
alteration in the electric potential in the cytoplasm, for both cells. 

Discussion 

There are other membrane phenomena due to electric charge 
effects that we have not considered in our model, such as the 
electric pump current [22]. We have not included these 
phenomena here as they seem to be negligible, in a first 
approximation, to the determination of the electric surface 
potentials. For example, Jackie [22] states that "the net electric 
pump current is not a relevant parameter in the casual theory of 
the resting potential". 

Experiments have revealed important discrepancies in the 
electrical properties of spinal ganglion neurons and neuroblastoma 
cells [9,16]. These differences are reflected in the parameter values 
we used in our simulations, as shown in Table (2). We notice that 
the values for the electrophoretic mobility, and the charge densities 
in the glycocalyx and on Sgg are much more negative for the 
neuroblastoma than for the ganglion. Because neuroblastoma cells 
contain a higher amount of negative pfg, Qs^.^, and compared 
with ganglion neurons, we indeed expected to obtain lower i^^^^ 
and (psi^^. values for the neuroblastoma in our simulations. 

Our simulations also demonstrate that variations in the electric 
charges frxed onto the inner surface of the membrane have a small 
effect on the electric potential of the surfaces that compose the 
neuronal membranes (Fig. (2)). We observe only a gende gradual 
drop m (j^Sj^^ of the spinal gangfion and the neuroblastoma cells 
during the resting state, as charges fixed on She decrease (more 
negative values). However, our model shows that decreasing the 
density of charges fixed onto proteins of the cytoplasm (increasing 
\pfi\ and Pfc/Pfg) results in a substantial decay of ijd^j , in both cells 
(Fig- (3)). 

Nonetheless, variations of intracellular charges fixed on the 
membrane and on the cytoplasmic proteins have no effect on the 
potential on Sgi, (Figs. (2) and (3)). This is related to the fact that 
the membrane plays a role in electrically isolating the intracellular 
and extracellular regions, due to the absence of charges within the 
lipidic bUayer (see Eq. (54)). 

The results we obtained for the spinal ganglion neurons and the 
neuroblastoma cells are generally similar to those obtained for the 
squid axon membrane in [4]. These authors showed that a 
decrease of 2si,. provokes a gentle decrease of i^^^^ of the squid 
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♦ — « 
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Figure 2. Sensitivity of the membrane surface potentials to inner surface charge density. Electric potential on the surfaces of regions of 
the membranes of the spinal ganglion neuron ( x ) and the neuroblastoma cell (*), as a function of the ratio QstJQs,^, as Qs,^ is kept constant. In the 
resting state (A), = — 190.97 mV, for the ganglion neuron and i^^^ = — 198.97 mV, for the neuroblastoma, when QstJQs,,^ = 1 (maximum values), 



while ( 



-194.25 mV (ganglion) and . 



-203.30 mV (neuroblastoma), for QstJQs.,=5^ (minimum). In the AP state (B), <l>s,^ =39.99 mV, for 



the ganglion neuron and ifig^ =29.99 mV, for the neuroblastoma, when Qsi^lQs,e = 1 (maximum), while il>s^^ = 'i9Sl mV (ganglion) and ^g^^ =29.97 
mV (neuroblastoma) for QstJQs,s~^^ (minimum). In all simulations, for resting and AP states, 1^5^^ = — 34.97 mV, for the ganglion, and 
^g^^ = —25.17 mV, for the neuroblastoma. In both graphs, Pf^ = 20pfg. 
doi:l 0.1 371 /journal.pone.00961 94.g002 



axon membrane. This behavior was also observed in our 
simulations for the membranes of spinal ganglion neurons and 
neuroblastoma cells. Nevertheless, their results indicate that a 
decrease of Qsi,^ causes a sensitive increase of (^^^^ during the AP 
state and a small decrease of during the resting state, whereas 
our results show a constant (l)g^^ value for ganglion and 
neuroblastoma cells. The insensitivity of (jig^^ to variations of 
Qstc which we have found seems more reasonable, given the above 
mentioned isolating effect of the Upidic bUayer. 

Cortez and collaborators [4] have shown that a decrease of py^ 
(in the same range of py^/ pj-g which we studied) causes practically 
no change in the surface potentials. A possible reason for this may 
be that the value for the squid axon is approximately zero, so 
that the values of py^, in the domain of their graphs, are very close 
to zero. In contrast, our simulations indicate that a decrease of py^ 
provokes an expressive fall of 05^^. In our case, (and py^) values 



for ganglion neurons and neuroblastoma cells are much more 
negative than those observed for squid axons and, therefore, a 
decrease of p^^, has a high influence on (j) . 

An interesting result of our calculations is that, in the spinal 
ganglion neuron, the electric potential across the glycocalyx 
decreases, and this does not occur in the neuroblastoma cell. This 
reveals an important discrepancy of the electric fields in the 
glycocalyx of both types of cells (Figs. (4) and (5)), and may explain 
the difference between their electrophoretic behavior, which was 
observed in experiments [9, 1 6] . As expected, the electric potential 
presents a Knear behavior within the bilayer of the membrane 
during the resting and AP states, due to the absence of electric 
charges in this region. 

The strong negative electric potential on S/,,. is a characteristic 
of the potential profile in the resting state, and this probably occurs 
for all types of neuronal cells (Fig. (4)). The steep increase of the 
potential from towards the bulk cytoplasmic region is regulated 



Action Potential 

— — Spinal Ganglion Neuron 
— • - Neuroblastoma 




0 10 20 30 40 50 

Pfc I PH 



0 10 20 30 40 50 



Figure 3. Sensitivity of the membrane surface potentials to charge density in the cytoplasm. Electric potentials i^^^^ and </i^^ as a function 

of Pfc/Pfg, as pfg is kept constant, for the spinal ganglion neuron ( x ) and the neuroblastoma cell (*). In the resting state (A), 4st^ = — 129.63 mV, for 
the ganglion neuron and ^g^^ = — 130.36 mV, for the neuroblastoma, when Qs^^lQs^ = 1 (maximum values), while ^^^^ = —213.74 mV (ganglion) and 
^Sj^ = —219.56 mV (neuroblastoma), for 2sk/2s,j =50 (minimum). In the AP state (B), (^^^^ = 82.26 mV, for the ganglion neuron and i^^^ =66.57 mV, 
for the neuroblastoma, when 2sj,/2s,j = l (maximum), while ^5^^=42.13 mV (ganglion) and 1^5^^=31.28 mV (neuroblastoma) for QstJQs,s = ^Q 
(minimum). In all simulations, for resting and AP states, (^^^ = —34.97 mV, for the ganglion, and <^s^ = —25.17 mV, for the neuroblastoma. In both 
graphs, Qs,, = 'i()Qs.,- 
doi:1 0.1 371/journal.pone.00961 94.g003 
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Figure 4. Electric potential across the membranes of spinal 
ganglion neurons and neuroblastoma cells, during resting 
state. Solutions of Eq. (52) with boundary , and Eq. (45) with 
boundary <j>+a3 = 'l>R result respectively in = —34.97 mV and 
= — 192.22 mV, for the spinal ganglion neuron (solid), and for the 
neuroblastoma cell (dashed) in ^ = —25.17 mV and i^^^^ = —200.66 mV. 
For all simulations, Qs^ = ^OQs^ and py^=20pj^. 
doi:l 0.1 371/journal.pone.00961 94.g004 

by the negative charges spatially distributed in the cytoplasm. 
Even though the net value of charges of proteins is predominantly 
negative in the cytoplasm, our simulations indicate that the 
contribution of these charges to the intracellular potential profile is 
much smaller than the effect of charges frxed on She- This is shown 
by the curvature of the potential in the cytoplasmic region. 

The neuroblastoma cells, like all cancerous cells, multiply 
quickly. Alterations of the dynamics of cellular multiplication 
mediate changes in the synthesis, structure and degradation of the 
membrane components [1 1], which result in deformations on the 
structure and composition of the surfaces of membranes [23]. 
These deformations provoke changes in the composition of electric 
charges on the membrane. Our results indicate that the alteration 
of these charges and of those within the cells may influence the 
behavior of the potential on the inner surface of the neuroblastoma 
cells. 

Experimental observations have suggested that the resting state 
and the generation of action potentials in human neuroblastoma 
cells depend on the degree of the morphological diflFerentiation of 
the cell. Some of these cells are relatively non-excitable [24,25]. 
Kuramoto et al. [26] stimulated the growth of SK-N-SH human 
neuroblastoma cells under standard culture conditions. These 
cancerous cells remained morphologically undifferentiated, par- 
tially responded to injections of pulses of electric current, and 
presented deficiency of the depolarizing component of the 
mechanism that generates the action potential. We included these 
findings in our simulations, and Fig. (5) is consistent with the fact 
that the depolarization of the electric potential in the neuroblas- 
toma, during generation of the action potential, is less intense than 
in the healthy spinal ganglion neuron. The neuroblastoma should 
thus generate a lower firing rate in response to its input excitation, 
and this may affect the transmission of signals through networks of 
these neurons and their functions of storage and communication of 
information. 

Mironov and Dolgaya [17] have suggested that the outer 
electric charges for the neuroblastoma cells and erythrocytes are 
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Figure 5. Electric potential across tKie membranes of spinal 
ganglion neurons and neuroblastoma cells, during AP state. 

Solutions of Eq. (52) with boundary (jd^ , and Eq. (45) with boundary 
il>+aD = 't'A result respectively in ,jig^ = -14.97 mV and (jig^^ =39.99 mV, 
for the spinal ganglion neuron (solid), and for the neuroblastoma cell 
(dashed) in ^i^^ = -25.17 mV and (l>s^ =29.99 mV. For all simulations, 
Qs^=iOQs,^ and pf,=20pfg. 
doi:1 0.1 371/journal.pone.00961 94.g005 

similar, but the spinal ganglion neurons strongly differ from these 
cells. Therefore, the molecular structure (and the resulting 
constitution of charges) on the outer surface of the membrane of 
the neuroblastoma cells would be similar to the erythrocytes, and 
may be constituted by ^ 40% of peripheral proteins and ^ 60% 
of gangliosides. Our results illustrate that the drop of the potential 
across the glycocalyx for the neuroblastoma cell is much smaller 
than for the spinal ganglion neuron, during both resting and AP 
states. This corroborates previous studies which show a smaller 
decay of the potential for the erythrocyte in the glycocalyx than for 
the neuron [2,4,5] . The different behavior of the potential across 
the glycocalyx, for the neuroblastoma and the spinal ganglion 
neuron, should indicate important differences among these cells, of 
the properties that enable the transmission of electric signals 
through the membrane. This occurs due to the fact that different 
molecular structures of these membranes interact differently with 
(i) the outer electric field, which is responsible for the orientation of 
the charged particles that are closer to the membrane, and (ii) the 
potential on the outer surface of the membrane. The nature of 
these interactions are crucial for many cell processes, such as the 
beginning of the process of triggering of the action potential, which 
depends on the opening of specific Na^ channels. 

Our results may also contribute to understanding the resistance 
of the neuroblastoma to certain chemotherapeutic treatments 
[27,28]. The smaller change of the potential, in response to 
changes in properties of cellular cultures (pH value's, for instance) 
and to the amount of fixed charges present in the membrane due 
to alterations in its composition and structure, may be an electric 
property responsible for the low pharmacological response. 
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